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Abstract. The eigenvalue equation of a band or a block tridiagonal matrix, 
the tight binding model for a crystal, a molecule, or a particle in a lattice with 
random potential or hopping amplitudes: these and other problems lead to 
three-term recursive relations for (multicomponent) amplitudes. Amplitudes 
n steps apart are linearly related by a transfer matrix, which is the product 
of n matrices. Its exponents describe the decay lengths of the amplitudes. A 
formula is obtained for the counting function of the exponents, based on a du- 
ality relation and the Argument Principle for the zeros of analytic functions. It 
involves the corner blocks of the inverse of the associated Hamiltonian matrix. 
As an illustration, numerical evaluations of the counting function of quasi ID 
Anderson model are shown. 



1. Introduction 



Consider the equation 



(1) 



where A^, B^, and C„ are given complex non singular square matrices of size m, 
i? is a parameter, and the vectors u„ S C™ are unknown. 

The case An = Aj^, C„ = b\_^ will be referred to as the Hermitian case, and is 
common occurrence in physics: it describes a chain of m— level atoms with Hamil- 
tonians An, and couplings Bn between neighbor atoms. In a lattice model for 
transport or in a model for a crystal, the atoms may be thought of as sections of 
the lattice or crystal, each containing m sites. In linear algebra, if off diagonal 
blocks are triangular, the equation represents the eigenvalue equation for a band 
matrix of bandwidth 2m -I- 1. 

At each n the equation provides u„+i in terms of u„ and u„-i, and an iterative 
evaluation of them can be started from initial conditions uq, ui. The process can 
be formulated via a 2m x 2m transfer matrix: 



(2) 



Un+l 
Un 



T{E) 



Ul 
Uq 



The n-step transfer matrix T{E) is the product tn{E)tn-i{E) ■ ■ ■ ti{E) of 1-step 
transfer matrices with the block structure (/„ is the unit matrix of size m) 



(3) 



tkiE) = 



B-\E-Ak) -B-'C, 
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We are concerned with the 2m exponents of the transfer matrix T{E), 

(4) £,a = -\0g\Zal 

n 

and their (normahzed) counting function: 

2m 
a=l 

An operative approach to the numerical evahiation of the exponents is found in 
ref.[T7j. In several cases of interest, the exponents have definite large n limits 
(Lyapunov exponents) and describe the exponential rates of growth or decay in 
n of the eigenvalues Za, and of the solutions m„ of the recursive eq.(IT]). They 
are important, especially the largest or the (positive) smallest, in the study of 
dynamical systems, stability, localization problems, wave transmission in layered 
structures |3j . Analytic expressions for the distribution of Lyapunov exponents are 
known only for few cases, and such cases do not originate from a Hamiltonian, i.e. 
a recursion equation like ([T]) [TOl [3 [2] . Most of the results where an Hamiltonian 
is given first, were obtained numerically [51 112[ [T5] or via some expansion [4,|lj. 
The main result of the paper is the following exact formula, valid for a single transfer 
matrix of type ([5]): 

Proposition 1.1 (Counting function). 

zG^(£;,z)i„B„- iG^(£;,z)„iCi , 
z 

z = e«+'"^ 

G^{E,z)ab are the blocks of size m of the resolvent G^{E,z) — [H^{z) — E]^^ , 
and H^{z) is the block tridiagonal matrix of size mn with corners 

Ai zBi Ci/z 
C2/Z -■■ -■■ 

zBn-l 

zBfi Cn/z An 

The suffix B stands for "balanced" , as this matrix is more tractable in numerical 
calculations than the Hamiltonian matrix H{z), to be introduced next. 
Notes: 1) the integral will be shown to be a contour integral, and cancels all z de- 
pendent terms of the Laurent expansion. Therefore one cannot remove the integral 
in the large n limit. 2) The imaginary part of the integral is zero. 
An expression of the counting function where the dependence on z is made "ex- 
plicit" will be given in the Appendix, prop. (15.11) . 



(7) 



H^iz) 



2. A SPECTRAL DUALITY 

According to the Argument principle of complex analyis, the zeros of an analytic 
function f{z) contained in the disk \z\ < R (inside the domain of analitycity) are 
enumerated by the following complex integral on a circle of radius R: 

dz f'{z) 
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2x2x80, w= 18 



Figure 1. The counting function A/'(^) (a different normalization 
is used) for the Anderson model in a cube 2 x 2 x 80, with disor- 
der strength w = 18, energy E = 0. Each discontinuity marks a 
positive exponent £,i, ■ ■ ■ ,£,4 (the negative exponents are just their 
opposites). 



We are interested in the zeros of f{z) = dct[T(£') — z], and their distribution Af{£,) = 
(l/2TO)n(e"^). The result for the counting function rests on a remarkable duality 
that relates f{z) to the characteristic polynomial of a Hamiltonian matrix H{z), 
which naturally arises from the problem ([l} restricted to n steps in a ring topology, 
but with generalized boundary conditions specified by a complex parameter z: 



(9) 



H{z) 



C2 ■■• 



zB„ 



Ci/z 



Bn-l 
An 



The polynomial f{z) — dei\T{E) — z] has degree 2m in the variable 2;, and degree 
nm in the parameter E. Duality states that the polynomial is proportional to the 
characteristic polynomial of the matrix H{z): 



Proposition 2.1 (Duality, P^Hi] ). 
(10) det[r(£;) - z] = 



, det[E-H{z)] 
det[Bi---BJ 



It is a duality relation as it exchanges the roles of the parameters z and E among 
matrices: whenever z is an eigenvalue of the transfer matrix evaluated at E , the 
value E is an eigenvalue of the block tridiagonal matrix H{z). 

Proposition 2.2. The matrix (z) is similar to H{z") : H{z'') = S{z)H^ {z)S{z)- 
where S(z) is the block diagonal matrix with blocks S{z)ab = SabZ°'Im- 
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As a consequence the two matrices have the same eigenvalues. Another consequence 
is the property {ze-'^''/") = S{e'^''/")H^ {z)S{e-'^''/"). 

Before deriving the expression for the counting function, we quote interesting and 
related analytic consequences of duality. Duality and Jensen's formula of complex 
analysis provide an equation for the exponents of the transfer matrix: 



Proposition 2.3 (13]) 

2m 

2m 



,=1 

(11) =— r"^log|det[i/(e<+'n-ii;]|-— ^logldetSfcl 

•'^ fc=i 

Corollary 2.4. For = 0, the formula gives the sum of the positive exponents: 

2m 

m ^-^ 

o=l 

(12) =— r''^log|det[77(e^'^)-i?]| Vlog|deti?,| 

fc=i 

Remark 2.5. The sum of all the exponents is i log | detr(£')|. Then: 

2m n 

(13) E = - E (i°g I I - log I I) 

a=\ k=\ 

Remark 2.6. // the transfer matrix is symplectic, i.e. there is a matrix E such 
that T{E)^YjT{E) — S, then the exponents come in opposite pairs. 

Proof. If z is an eigenvalue of T{E) with eigenvector 9, then S0 is an eigenvector 
of T{E)'^ with eigenvalue 1/z . This means that 1/z is eigenvalue of T{E). The 
corresponding exponents are either both zero or opposite real numbers. □ 



Eq. (jl2[) is exact and applies to any single transfer matrix. However it is rem- 
iniscent of the formula for the sum of the Lyapunov exponents of an ensemble 
of random transfer matrices, obtained by Herbert, Jones and Thouless for 2 x 2 
matrices, and by Kunz, Souillard, Lacroix |llj for larger ones: 

— y^\a{E)= I dE'p(E') log \E - E'\ + const. 

771 -^-^ I 

This formula was obtained for the Hermitian case, but holds also in non Hermitian 
tridiagonal problems [SI [9] . In eq. ([T2l) the angular average on replaces the en- 
semble average, that produces the ensemble density of eigenvalues. It seems that 
an ergodic property is at play for such systems. 



3. The counting function 



The expression ([7]) for the counting function can be obtained straightforwardly 
from eq. (fTT|) . Here we give a direct proof based on the famous formula eq.®, the 
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duality relation, and the simple formula 



4- det M{z) = det M{z) tr 
dz 



M{z)-^-^M{z) 
dz 



Proof. Let G{E,z) = [H{z) - E]-^, f{z) = det[r(£;) - z], and use duality: 



£[z-^det[H{z)^E]] _m ^ 



z"'det[H{z)-E] 



G{E,z)-H{z) 



m 



+ tr 



G{E, z)inBn nG{E, z)nlC\ 



The complex integral of the first term is m. Inclusion of the normahzation factor 
gives 



dz 



2 2m /|2|=en5 2-Kiz 



tr 



Z G{E, z)xnBn G{E, z)niCi 



1 1 

2 2m 



d0 



-tr 



Z G{E, z)inBn G{E, z)nlCi 

z 



where z = e"?+*^. Next wc set 9 = nip, then z = e"'^^+*'^\ We introduce the 
resolvent of H^{z) and note that G{E,w"')ab = w°-~''G^ {E,w)ab by virtue of the 
similarity; in particular it is: 



w"G{E,w")in = wG''{E,w)i„, 



■\G{E,w^)ni = -G''{E,wU 



The final formula is obtained (after renaming k; as ^;). Since the counting function 
is real, only the real part of the integral is nonzero, and: 



(14) = / " ^ Im tr 



z G^(£, z)inBn --G"" (E, zUCi 
z 



□ 



Remark 3.1. While the Argument Principle requires the circle radius to be e"^, 
the use of the balanced martix scales the radius to e^ and makes the formula for the 
counting function useful for numerical studies. 

In numerical calculations the angular integral is replaced by a sampling on dif- 
ferent angles. As the radius crosses an eigenvalue e^" the counting function 
jumps to a new plateau after a strong oscillations due to close sampling points: 

The formula requires the inversion of the matrix , and the explicit depen- 
dence on the interesting parameter z is lost. In the Appendix wc show that this 
dependence can be made more explicit by means of the resolvent equation for the 
matrix with corners removed. 



6 LUCA GUIDO MOLINARI AND GIUSEPPE LACAGNINA 



3x3x80 




0,5 1 1.5 

5 



Figure 2. The counting function A/'(^) (a different normalization 
is used) for the Anderson model in a cube 3 x 3 x 80, with disor- 
der strength w = 18, energy E = 0. Each discontinuity marks a 
positive exponent ^i, . . . , ^g. 

4. Anderson Model 

As an illustration of the formula, consider the Anderson model for the motion of 
a particle in a cubic lattice with the shape of a bar, with a random potential. The 
particle may hop with equal amplitudes from a site k to neighboring ones k + e, 
and is subject to a site-potential e(fc) (a random number uniformly distributed in 
the interval [—w/2,w/2]): 

'^ij{k + e} + e{k)%lj{k) = Ei}{k) 

e 

If the bar is sliced perpendicularly to the long direction [k — k^), each slice is a 
square 2D lattice, with m sites. Accordingly the Hamiltonian gains a block struc- 
ture, with diagonal blocks describing slices, and blocks Bk = Ck = Im describing 
hopping among neighboring slices. The Schrodinger equation is now written as in 
([T]): Uk+i + Uk-i + AkUk — E Uk, where the m components of Uk are the amplitudes 
ip{k) in the slice k. 

Because of disorder, in the large n limit, eigenstates are exponentially localized 
(Anderson localization) on a distance given by the inverse of the smallest positive 
exponent. In the large n limit all the exponents ±^a(i?) exist, and are independent 
of the realization of the disorder. 

The two figures show the counting function for a single realization of disorder, and 
a sample of length 80 (cross section 2 x 2, or 3 x 3). Note the decrease of distances 
among exponents. 

The subject on Anderson localization is vast. A discussion of the spectral prop- 
erties of Anderson matrices with non Hermitian boundary conditions (the z and 
l/z dependence) is presented in [IS]. The one dimensional case was studied by 
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Goldsheid and Khoruzhenko (Hermitian case) [8j and several others, and by us in 
the non Hermitian case 1161. 



(15) 



5. Appendix 

Let h be the matrix H(z"') with corners removed (then it is z-independent): 

C2 ■■• ■■• 



h = 



and let g{E) = [h — E] ^ be the resovent matrix, with blocks gab- Then: 
Proposition 5.1. 



(16) 



tr 



zG^(£;,z)i„B„--G^(S,z)„iCi 



tr 



— Cignn —C\gnl + In 



Proof. A(z") — H{z") — /i is a matrix whose non zero blocks are A(z")i„ = z "Ci 
and A(z")„i = z'^Bn- The matrix equation for the resolvent is 

(17) G{E, z") = g{E) - g{E)^{z^^)G{E, z") 

The four equations for the blocks a,b — l,n can be put in matrix form 



Gin Gil _ gin gii _ gin gii ^"S„ 

Gnn Gnl g^in gnl gnn gnl 

Gi/z" 

The blocks of G{E, z") are expressed in terms of the blocks of G^{E, z) 



Gin 
Gnr 



Gil 
Gnl 



Gin Gil 
Gnn Gnl 



1/z" 
z" 



zGfn ^"Gfi 



Then 



Z'^ln Z 
^nn/^^^ Gnll'^ 

z'^gin 



Z"5i„ Z"gii 
gnn/z"^ .9nl/z" 



z^'gii 




' Bn ■ 




gm/z'^ 







Cl . 






-1 




9in gii 




■ 1/z" 


Qnn Qnl 








Z^ln Z (jii 
^nn/z^ Gnl/z 





' Bn ■ 




+ 


Ci 





with solution 

„riB ./n/^B 
Z'-'ln Z i^ii 

^nn/z"' Gni/z 

Finally, left multiplication by the diagonal matrix with blocks _B„ and — Ci gives 

zBnGin Z^ BnGii 

-GiGnn/z^ —CiGni/z 

-1 



9ln 911 
Qnn Qnl 



B-^/z- 






^z"cr 



Im 

-/„ 
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The r.h.s. is simplified by means of tfie following identity for invertible matrices A 
and B: {A-^ + B"!)-! =A-A{A + B^^A, 



' Irn 




' Im " 


. -Ira 




-Im _ 



^"Bn 




gin 


911 




' Im ■ 


r 


' Im " 


-^Ci . 




9nn 


9nl 


+ 


. -Im _ 




. -Im _ 



The trace is then taken: 



(18) tr 



zBnGf^ CiG^i 
z 



tr 



-z''Bn9ln - Im. -^"-Bn5ll 
^Cignn prCl5nl + In 



As a further step one could perform the the block inversion of the matrix by Schur's 
formula. Since the trace is then taken, only the diagonal blocks are needed: 

-1 



= -tr 



z^Bngin + Im + -B„fifll(— Ci5„i + Im) ^Cig„ 

2" 



(19) 



+ tr 



— Cignl + Im + Cignniz" Bngin + Im) ^-Bnffll 



□ 



This final expression allows in principle the evaluation of the counting function 
with a single big matrix inversion for obtaining the four corner matrix blocks gab 
{a,b = l,n). Then, for each ^ value, the inversion of a square m x m matrix is 
required. 



References 

[1] H. Schulz-Baldes, Perturbation theory for Lyapunov exponents of an Anderson model on a 

strip, Gcom. and Funct. Anal. 14 (2004) 1089-1117. 
[2] C. W. J. Bccnakkcr, Random matrix theory of quantum transport, Rev. Mod. Phys. 69 

(1997) 731-808. 

[3] A. Crisanti, G. Paladin and A. Vulpiani, Products of Random Matrices in Statistical Physics, 
Springer, Berlin (1993). 

[4] B. Derrida, K. Mecheri and J. L. Richard, Lyapunov exponents of products of random 
matrices: weak disorder expansion. Application to localization, J. Physique 48 (1987) 733- 
740. 

[5] B. Derrida, J. L. Jacobsen and R. Zeitak, Lyapunov exponent and density of states on a 

one- dimensional non-Hermitian Schrodinger aquation, J. Stat. Phys. 98 (2000) 31-55. 
[6] T. Kottos, A. Politi, F.M. Izrailcv, and S. Ruff'o, Scaling properties of Lyapunov spectra for 

the band random matrix model, Phys. Rev. E 53 (1996) R5553. 
[7] A. Gamba and I. V. Kolokolov, The Lyapunov spectrum of a continuous product of random 

matrices, J. Stat. Phys. 85 (1996) 489. 
[8] I. Ya. Goldsheid and B. Khoruzhenko, Eigenvalue curves of asymmetric tridiagonal random 

matrices, Electronic Journal of Probability 5 (2000) 1-28. 
[9] I. Ya. Goldsheid and B. Khoruzhenko, Thouless formula for random non-Hermitian Jacobi 

matrices, Isr. J. of Math. 148 (2005) 331-46. 
[10] M. Isopi and C. M. Newman, The triangle law for Lyapunov exponents of large random 

matrices. Comm. Math. Phys. 143 (1992) 591-598. 
[11] J. Lacroix, Computations of the sum of positive Lyapunov exponents for the Lloyd model in 

a strip, in Lyapunov Exponents, Lect. Notes in Mathematics 1186, Springer- Verlag (1984). 
[12] P. Markos and B. Kramer, Staistical properties of Lyapunov exponents and of quantum 

conductance of random systems in the regime of hopping transport, Ann. Physik 2 (1993) 

339-360. 



COUNTING THE EXPONENTS OF SINGLE TRANSFER MATRICES 



9 



[13] L. Molinari, Spectral duality and distribution of exponents for transfer matrices of block 
tridiagonal Hamiltonians, J. Phys. A; Math. Gen. 36 (2003) 4081. 

[14] L. G. Molinari, Determinants of block-tridiagonal matrices, Linear Algebra and its Appli- 
cations 429 (2008) 2221-2226. 

[15] L. G. Molinari, Nan Hermitian spectra and Anderson Localization, J. Phys. A: Math. Theor. 
42 (2009) 265204 (15 pp). 

[16] L. G. Molinari and G. Lacagnina, Disk-annulus transition and localization in random non- 
Hermitian tridiagonal matrices, J. Phys. A: Math. Theor. 42 (2009) 395204 (9 pp). 

[17] K. Slevin, Y. Asada and L. I. Deych, Fluctuations of the Lyapunov exponent in two- 
dimensional disordered systems, Phys. Rev. B 70 (2004) 054201 (7 pp). 

[18] Y. Y. Zhang and S. J. Xiong, Statistics of Lyapunov exponents of quasi one dimensional 
disordered systems, Phys. Rev. B 72 (2005) 132202 (4 pp) 



Physics Department, Universita degli Studi di Milano and I. N.F.N, sez. Milano, Via 
Celoria 16, 20133 Milano, Italy 

E-mail address: luca . molinariSmi . inf n . it 



